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ABSTRACT 

The velocity dispersio n of stars in the solar ne ighbourhood thin disc increases with 
time after star formation. INordstrom et al.l {2004) is the most recent observational at- 
tempt to constrain the age-velocity dispersion relation. They fitted the age-velocity 
dispersion relations of each Galactic cardinal direction space velocity component, U 
(towards the Galactic centre), V (in the direction of Galactic rotation) and W (to- 
wards the North Galactic Pole), with power laws and interpreted these as evidence 
for continuous heating of the disc in al l directions througho ut its lifetime. We re-visit 
these relations with their data and use iFamaev et al.l (|2005f ) to show that structure in 
the local velocity distribution function distorts the in-plane (U and V) velocity distri- 
butions away from Gaussian so that a dispersion is not an adequate parametrization of 
their functions. The age-aw relation can however be constrained because the sample 
is well phase-mixed vertically. We do not find any local signature of the stellar warp 
in the Galactic disc. Vertical disc heating does not saturate at an early stage. Our 
new result is that a power law is not required by the data: disc heating models that 
saturate after ~ 4.5 Gyr are equally consistent with observations. 
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The velocity dispersion of the molecular clouds in which 
star formation is taking place is very much less, by a factor 
of three to five, than is the (one-dimensional) dispersion of 
old stars near the Sun. The (set of) processes which change 
the velocities in any direction of initially slow-moving young 
stars into 'hotter', higher-dispersion orbits with time are disc 
'heating' mechanisms and are evidently fundamental aspects 
of disc galaxy evolution. 

The most recent observational atte mpt to constrain 
the ag e-velocity dispersion relation is by INordstrom et all 
(2004), hereafter N04, using their Geneva-Copenhagen sur- 
vey of the solar neighbourhood. This is the largest homoge- 
neous data set available with space velocity components, U 
(towards the Galactic centre), V (in the direction of Galac- 
tic rotation) and W (towards the North Galactic Pole), and 
isochrone-derived stellar ages. N04 fitted the age-velocity 
dispersion relations of each space velocity component (and 
the total) with continuous smooth power laws: a oc age fe , 
where the scaling (heating) exponent, k, is 0.31, 0.34, 0.47 
and 0.34 for ajj , ay, o~w and a to tai respectively. They inter- 
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preted these relations as evidence for continuous heating of 
the disc in all directions throughout its lifetime. 

However, continuous vertical (aw) heating is unex- 
pected. Stars spend most of their time near the apocentre of 
their vertical orbits, furthest away from the Galactic plane. 
For stars to experience continuous vertical heating through- 
out their lifetimes requires heating processes that operate 
away from the plane. G iant molecular clouds (GMCs) gravi- 
tationally scatter stars ( Spitzer fc S chwarzschild 1951, 19531 : 



lLacevlll984l : IVillumsenl 



19851 : lHanninen fc Flvnnl 120021 ) but 



they exist in the extreme thin disc. This suggests that ver- 
tical heating should not continue throughout the lifetime of 
the disc but should saturate after some epoch if GMCs are 
indeed the primary cause. 

There are many heating mechanisms in the literature 
using non-disc pr ocesses to explain contin uous vertical heat- 
ing observations. lLacev fc Ostrikerl l|l985h proposed scatter- 
ing by m assive black holes in the halo, which have been sim - 
ulated bvlFuchs et al.l l|l994h and lHanninen fc Flvnnl (|2002l ). 
ICarr fc Lacevl (|l987f l proposed dark clusters of less massive 
objects could heat stellar discs. However, neither of these 
major mechanisms have any direct observational support. 
Minor mergers with dwarf galaxies is another possible heat- 
ing mechanism dToth fc Ostrikerl Il992l : IWalker et a"illl996l : 
iHuang fc Carlberd ll997F T" Dramatic disc heating can oc- 
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cur when massive satellite galaxies fall on to galactic discs 
( Oumn et al.ll 19931 ). which could leave a sharp feature in the 
a,ge-aw relation but many smaller galaxies acc reting could 
heat a disc smoothly l| Velazquez fc Whitdl 19991 ) . Other pos- 
sible vertical heating mechanisms that may play a minor 
role are the outer Lindblad resonance from the Galactic bar 
l|Kalnaislll99TI ; |Pehnedfl999l l2000l;lFuxll200ll ; lOuilledliool ) 
and evaporating star clusters (|Kroupal l2002). In this paper, 
we re-visit the age-o^ relation using the same data as N04 
to test the uniqueness of their power law model fit. We in- 
vestigate whether the data permit models where the vertical 
disc heating saturates, thus negating the need to invoke hy- 
pothetical or poorly understood components of the Galaxy. 

Whether heating is continuous or saturates with 
time requires separate consideration in the plane (au 
and ay) and o u t of the plane (aw)- Simulations by 
|Pe Simone et all |2004) illustrate that transient spiral 
waves are most likely to be the dominant in-plane heat- 
ing mechanism in the solar neighbourhood. Transient spi- 
ral waves can lead to gravitational potential fluctua- 
tions in the disc that excite random motions of disc 
stars dBarbanis fc Wolfe Il967l; ISellwood fc Carlberdjl984l ; 



ICarlberg fc Sellwoodl Il985l ~l Sellwood fc Binnevl 12002) but 
only effectively in U and V because they are confined to 
the thin disc. GMCs also contribute to the in-plane heat- 
ing but are unlikely t o be the dominant heating source 
l|De Simone et aLll2004l ). Once a star's epicyclic amplitude 
reaches a size-scale similar to that of the spiral arm radial 
wavenumbers, the star will cross more than one spiral arm 
as it radially oscillates m aking it more difficult for th e spi- 
ral structure to perturb it |Carlberg fc Sellwood| [l985). This 
means spiral heating is theoretically expected to become 
much less efficient with time so that its in-plane heating 
effectively saturates. 

However, whether au and ay should actually saturate 
requires simulating the heating from both realistic Galactic 
spiral structure and GMCs but these calculations are com- 
plicated. Spiral features heat a disc at a rate that depends on 
the spatial pattern of the potential and on its time variabil- 
ity. In principle, the photometry of real galaxies can yield the 
spatial structures of galaxies and temporal information con- 
cerning spiral structure can be determined from numerical 
simulation and dynamical theory. In practice, the strengths, 
growth rates and duty cycles of spirals are not well con- 
strained. Therefore, the relative effectiveness of spiral and 
GMC heating has to be modelled with an e mpirical param- 
eter i|jenkins fc Binnevlll990l ; I Jenkins! H992T ). Naturally, the 
empirical parameter was guided by observation and so the 
models reflected the only observationa lly determine d age- 
velocity dispersion relation of the time l|Wielenll 19771) . 

Like N04, IWielenl (|l977l ) also found a continuous power 
law for a total (and aw) with k ~ 0.5 but this was based on 
ranking the age of mainly K and M dwarfs using chromo- 
spheric emission l ine measurements, wh ich is fraught with 
systematic errors. ICarlberg et ail |l985l) used more reliable 
isochrone ages and found observational evidence for satu- 
ration in the age-a tota i relation at ages > 5 Gyr but this 
was obtained using tange ntial velocities and only 255 stars. 
lOuillen fc Garnettl (|200ll ) argued for earlier saturation in 
au, ay, aw and a tota i between 3 and 10 Gyr. H owever, this 
was based on only 189 lEdvardsson et al.l (| 19931 ) stars (~ 21 
stars per logarithmic age bin) and although they used new 



improved isochrone ages from iNg fc Berteliil lll99Sl). their 
age de t ermination is les s reli able than N04 . Also, Stromgrer] 
i| 19871 ) . lFreemanl (|l99ll ) and lGomez et all i| 19971 ) all argued 
for a turnover in the observational relation, indicative of 
saturation. 

The in-plane heating saturation issue is complex be- 
cause it is not clear if the observations are reliable and the 
simulations depend on the observations. In this paper, we ar- 
gue that this issue is still not observationally resolved even 
though N04 presented evidence of continuous in-plane heat- 
ing. We show their analysis is not conclusive because the 
interpretation of the in-plane age-velocity dispersion rela- 
tion is complicated by the presence of strong, stellar age 
dependent, small-scale structure. 

|Pe Simone et all (|2004l ) demonstrate that heating by 
strong transient spiral waves can induce strong small-scale 
structure in the local U and V velocity distributions by im- 
parting momentum to stars to put them in specific regions 
of the [T V-plane, creating streams w ith a stellar age range of 
~ 3 Gvr. lSellwood fc Binnevl l|2002l ) found that spiral waves 
cause radial migration in the galactic disc near their coro- 
tation radius, preserving overall angular momentum leading 
to only a small increase in random motions. This naturally 
explains the 'dynamical stre ams' found in the kin ematic sur- 
vey of local K-M giants by iFamaev et all (|2005l ). hereafter 
F05, which have similar properties to moving groups and su- 
perclusters but which contain stars spanning a much wider 
age range. It suggests that the F05 streams are due to an 
inhomogeneous Galactic potential perturbing smooth star 
formation rather than the streams resulting from inhomo- 
geneo us star formation in a smooth potential. IFamaev et alj 
i|2007l ) show that both situations exist in the N04 sample: 
the Hyades-Pleiades stream (see Fig. |3} is mainly composed 
of field-like stars but also partly of coeval stars evaporated 
from the primordial Hyades cluster (~ 85 and 15 per cent 
of the stream for low-mass stars respectively). Alternatively, 
lOuillen fc Minchevl l|2005h propose the Hyades-Pleiades and 
Coma Berenices (or local) streams (middle branch in the 
top left panel of Fig. |3J) correspond to nearly closed orbits 
trapped at the 4:1 inner Lindblad resonance of a two-armed 
spiral density wave. 

IPehnenl (1 19981 ) used a maximum likelihood algorithm 
to estimate the velocity distribution of the N04 stars before 
their radial velocities were available using only parallaxes 
and tangential velocities. He found the traditionally recog- 
nised moving groups as well as a smooth background distri- 
bution, akin to the Schwarzschild ellipsoidal model, as did 
F05. The aforementioned streams and the Sirius dynami- 
cal stream all share the same region of UV phase-space as 
the Galactic thin disc background (see Fig. [3j). The result- 
ing non-Gaussian U and V velocity distributions pose the 
question: what parametrization of a complex distribution 
function is adequate to represent a 'dispersion'? 

The Gaussianity of the age-aw relation permits it to be 
re-visited in this paper. We investigate its derivation from 
N04 data by checking the stellar warp of the Galactic disc 
does not contribute to aw and that W is sufficiently well 
mixed that the presence of streams does not cause aw of 
the whole sample to deviate radically from its value for 
the smooth background (as is certainly happening in au 
and ay). This is demonstrated by showing that dynamical 
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streams do not affect aw by deriving the a,ge-aw relation 
with and without the Hercules stream. 

The Hercules stream is the only stream in the N04 sam- 
ple that does not share the same region of UV phase-space 
as the Galactic thin disc background, allowing it to be un- 
ambiguously defined in UV phase-space and removed in W 
without also excluding a large region of the background ellip- 
soid (as would happen in removing one of the other streams) . 
The Hercules stream has been speculated to be a signature of 
the Galactic bar. Whether Hercules stream stars are chaoti- 
cally, gravitational ly sca ttered from the inner Galaxy by the 
Galactic bar (|Fuxll200lh or they are more local stars cou - 
pled to the outer Lindblad res onance (|Dehnenlll999l . l200Ch 
is still uncertain. iQuillenl {2003) confirms that the Hercules 
stream remains a strong feature of the local velocity distri- 
bution function whe n the effect of s piral s tructure is added 
to the Galactic bar. iBensbv et al.l (|2007t ) investigated the 
abundance and age structure of the Hercules stream. They 
found it contains a mixture of young and old stars, all with 
UV kinematics more characteristic of the thick disc than the 
thin disc. 

Section[2]describes the N04 data used in this paper: sur- 
vey design (Section 12. space velocities (Section 12. 2[) and 
discusses the N04 age derivation (Section I2.3p . Our analy- 
sis is presented in Section where we reproduce the N04 
age-velocity dispersion relation age bins (Section I3.1|l and 
physically interpret them by qualitatively comparing them 
to the F05 velocity distributions ( Sections 13 . 21 and 13 . 3p . Sec- 
tion [3]4] confirms the stellar warp of the Galactic disc does 
not contribute to aw- We re- visit the age-ow relation us- 
ing the same age binning as N04 in Section 13.51 and with 
our own higher resolution age binning in Section \3. 61 before 
discussing our model fitting results in Section 2] 



2 DATA 

2.1 Survey design 

N04 presented a complete, all-sky, magnitude-limited, and 
kinematically unbiased sample of 16,682 nearby F and G 
dwarf stars. F and G dwarf stars are good probes of Galactic 
evolution because they are numerous and long-lived. Their 
radial velocities can be measured accurately and their ages 
can be estimated by comparison with stellar evolution mod- 
els, at least for the more evolved stars. Ages are crucial to 
place the kinematic properties of stars in an evolutionary 
context. The following section summarizes the properties of 
N04 most pertinent to re- visiting the age-velocity dispersion 
relation. 

The N04 observational input catalogue was selected 
from a compilation of catalogues available in the literature 
with Stromgren uvby/3 pho tometry of nearby F and G stars, 
mainly from the surveys bv lOlsenl l|l983l . ll993l . ll994al lbl) . The 
final sample is volume complete to « 40 pc. It is magnitude 
complete to V ^ 7.7 for the bluest stars (slightly fainter for 
the reddest G stars) and has a V cut-off magnitude of m 8.9 
(and 9.9 for the reddest G stars). 



2.2 Space velocities 

The input catalogue was observed using the CORrela- 
tion RAdial VELocities (C ORAVEL) phot o electr ic cross- 
corre lation spectrometers l|Baranne et al.l Il979l ; iMavorl 
I1985T I . Two CORAVELs cover the entire sky between them: 
one is on the Swiss 1-m telescope at Observatoire de 
Haute-Provence, France, and the other on the Danish 1.5- 
m telescope at ESO, La Silla. Their fixed, late- type cross- 
correlation template spectra match the spectra of the ma- 
jority of the input catalogue stars. The multi-epoch radial 
velocities (generally two or more spectroscopic observations) 
have a typical mean error of 0.5 km s _1 or less. 

The vast majority of t he N04 stars hav e proper motions 
in the Tycho-2 catalogue |H0g et al . 2000). This catalogue 
was constructed by combining the Tycho star-mapper mea- 
surements of the Hipparcos satellite with the Astrographic 
Catalogue based on measurements in the Carte du Ciel and 
other ground-based catalogues. A small number of stars have 
only one measurement, either from Hipparcos or Tycho. The 
typical mean error in the total proper motion vector is 1.8 
mas yr _1 . 

The primary source of dista nce fo r N04 stars is Hip- 
parcos trigonometric parallax (ty, ESA 1997), from which 
absolute magnitude (My) can be derived. This is adopted if 
its relative error (a^/n) is accurate to 13 per cent or better, 
otherwise th e photome t ric di stance calibrations for F and 
G dwarfs by ICrawfordl (| 19751 ) and lOlsenl (|l984l ) are used, 
with an uncertainty of only 13 per cent. Distances are not 
provided for stars with unreliable Hipparcos parallaxes (> 
13 per cent) and when photometric distances cannot be cal- 
culated. This occurs when the star is missing the necessary 
photometry and/or it falls outside the photometric calibra- 
tions. The absence of a distance estimate or radial velocity 
measurement reduces the size of the N04 catalogue with all 
full six dimensional phase-space information to 13,240 stars. 
The space velocity components (U, V, W) are computed for 
all the stars with (mean) radial velocities, proper motions 
and distances. Sp ace velocities in the Viz ieR catalogue ac- 
companying N04 l|Nordstrom et all l2004h are provided to 
the nearest km s _1 . The average error in each component is 
1.5 km s . 



2.3 Ages 

2.3.1 N04 age derivation 

N04 derive isochrone ages for their stars. Ages are derived 
by using the photometrically measured effective tempera- 
tures (T c ft), V magnitudes and metallicity ([Fe/H]) to place 
stars in the Hertzsprung-Russell diagram and reading off 
the age of the stars by interpolation between theoretically 
computed isochrones. T c s, V and [Fe/H] all require extinc- 
tion correction. The Stromgren /3 photometry provides in- 
terstellar reddening, E(b — y), fo r F stars via the intrinsic 
colour calibration bv lOlser] ( 1988h . This has been applied in 
the photometric temperature and distance determination if 
E(b — y) 0.02 and the distance is above 40 pc, otherwise 
the stars are assumed to be unreddened. T e g is calculated 



from the reddeni ng-corrected b — y, ci and mi indices and 
the calibration of lAlonso et al.l l|l996h . which is based on the 
infrared flux method. 

The majority of N04 metallicity values are derived 
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using the F or G star calibrations of ISchuster fc Nissenl 
l| 19891 ). Due to the ve ry few spectroscopi c calib rators avail- 
able at the time, the ISchuster fc Nissenl (|l989l ) calibration 
yields systematic errors for the reddest G and K dwarfs 
(b — y > 0.46). To remedy this, N04 performed a new fit 
of the uvby indices to [Fe/H] derived from high-resolution 
spectra for 7 2 dwarfs (0.44 ^ b — y ^ .59), using the same 
terms as the ISchuster fc Nissenl (Il989|) G star ca l ibrati on. 
When valid, N04 adopted the [Edvardsson et"al] (I1993T ) /3 
and mi calibrat ion for stars with high T e g and low gravi- 
ties outside the [Schuster fc Nissenl (Il989l ) calibrations. For 
stars outside both calibrations, N04 used other spectroscopic 
sou rces to derive a new rela tion using the same terms as 
the ISchusterfc^ Nissenl (|l989l ) F star calibration. The typi- 
cal [Fe/H] uncertainty is ~ 0.1 dex. 

Individual stellar ages are determined by a Bayesian es- 
timation method. For ev ery point in a dense grid of interpo- 
lated Padova isochrones (|Girardi et al.| [2000; Salasnic h et al.l 
2000) and given a star's nominal position in the three dimen- 
sional Hertzsprung- Russell 'cube' defined by log T e ft, Mv 
and [Fe/H], N04 compute the probability that the star is 
actually located at every po int. This procedure takes int o 
account observational errors |j0rgensen fc L indcgrcn 2005). 
Integrating over all the grid points gives the global likelihood 
distribution for the possible ages of the star. N04 account for 
statistical biases in the integration. The resulting probability 
distribution function ('G-function') is normalised to unity at 
maximum (see their fig. 13). A stable age estimate without 
significant bias is given by the maximum of the smoothed 
G-function. Upper (o-^g| h ) and lower (0^°™) 1<t confidence 
limits are the points where the G-function reaches a value 
of 0.6. 

Isochrone ages can only be determined in practice for 
stars that have evolved significantly away from the Zero Age 
Main Sequence. Therefore, unevolved (i.e. relatively young) 
stars in the N04 catalogue do not have reliable isochrone 
ages. This does not affect this study as we are concerned with 
the functional form of the age- velocity dispersion relation at 
the evolved end (i.e. old stars) and do not require a complete 
or unbiased age distribution function. 



2.3.2 Age comparisons between N04 and other data sets 

Holmbcrg et al. I (|2006h substantially improved temperature 
and metallicity calibrations of the photometry used by N04, 
evaluated their effects on the computed ages and performed 
extensive numerical simulations to verify the robustness of 
the derived relations. They found good agreement between 
their new ages and those of N0 4, supporting the N04 results. 

However. lHavwo"od" d2002f) demon s trates that the orig- 
inal calibration by ISchuster fc: Nissenl (|l989l ) is subject to 
systematic errors, mainly due to incorrect placement of the 
Hyades reference s equen c e. Using t h e sam e formalism as 
ISchuster fc Nissenl (| 19891 ). iHavwoodl (|2002h derived an al- 
ternative calibratio n, which is in good agreement with spec- 
troscopic data sets. Reid et al. (120071 ) calculated the [Fe/H] 
difference between IHavwoodl ( 2002) and N04 and found sys- 
tematic differences b etween -0.0 8 1 and -0.015 dex for dif- 
ferent colour ranges. iReid et all (|2007t ) al so compared the 
photometrically-derived [Fe/H] of N04 and IHavwoodl (|2002l ) 
with high-resolution echelle spectroscopically-derived [Fe/H] 
obtained by the Keck, Lick and AAT planet search pro- 



grams (IValenti fc Fischer! 120051 ). They found that while the 
dispe rsions were simila r (~ 0.1 dex), the differe nce be- 
tween Iwjentnnnichei] l|2005l ) and lHavwoodl (I2OO2I) (-0.041 
to +0 .045 dex) is smaller than between IValenti fc Fischerl 
|2005l ) and N04 (0.015 to 0.090 dex). 

Since the N04 ages are derived from isochrone fitting, 
systematic [Fe/H] errors translate to systematic age offsets. 
As a recalculation of the N04 ages to investigate the effect 
of systematic [Fe/H] errors is unnecessary for this paper, we 
proceed with the assumption that any age offset will have a 
negligible impact on our statistical use of age in this study. 
In addition, we are only using single stars with well-defined 
ages and age errors, so we assume the N04 age errors are the 
appropriate errors. 



3 ANALYSIS 

3.1 Reproducing the N04 age- velocity dispersion 
relation age bins 

N04 caution that many stars in their sample are binary or 
multiple systems, for which the derived metallicities and 
thus ages will be unreliable. They provide a catalogue flag 
(fb) that identifies confirmed or suspected binaries, where 
the information can come from one or several sources such 
as photometry, radial velocity or astrometry. N04 used only 
single stars (as defined by a null ft catalogue entry, 8,589 
stars) with relative age errors < 25 per cent (2,852 stars) 
to derive their age-velocity disperson relations in their fig. 
31. Our attempt to reproduce this sample by rejecting stars 
with fb catalogue entries to define single stars and selecting 
stars with (cr^gf 11 - age)/age < 0.25 and (age - oig^)/age ^ 
0.25 yields a smaller sample (2,801 stars). The missing 51 
stars are presumably due to N04 using more precise ages 
than their catalogue ages, which have been rounded up to 
the nearest 0.1 Gyr. 

N04 split their sample of single stars with relative age 
errors < 25 per cent into ten age bins of equal numbers. 
Due to the presumably lower age resolution in their cata- 
logue, we are unable to make the numbers of stars in each of 
our age bins exactly the same. This means that the bound- 
aries of our age bins are only similar to those used by N04, 
rather than being exactly the same. Unlike N04, we plot the 
space velocity diagrams (Fig. [TJ and histograms (Fig. |2J) of 
a representative age bin (plots of all the age bins are in the 
Appendix). The 3 km s -1 velocity bin sizes in Fig. [2] (and 
Fig. IA2[I are chosen to be twice the average error in each 
component of the N04 space velocities so that substructure 
can be resolved (F05 found the Hyades-Pleiades and Sirius 



dynamical streams have ov ~ 5 km i 



From the apparent 



lack of substructure in Fig. [T] (and Fig. IA1|I , it is not clear in 
Fig. [2] (and Fig. lA2[) whether the peaks in the histograms are 
Poisson noise, due to the relatively small number of stars in 
each bin at this velocity resolution, or substructure. Decom- 
position of the N04 sample into its constituent kinematic 
groups, like F05, is beyond the scope of this paper. Instead, 
we qualitatively compare the N04 and F05 space velocity 
diagrams and histograms to assess whether we expect com- 
plex substructure in all the age bins in Figs. IA1I and IA2I 
or whether the peaks are noise such that the distribution 
functions can be approximately described by Gaussians. 
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Figure 1. U — V (left), U — W (middle) and V — W (right) space velocity diagrams of the 278 N04 single stars with relative age 
errors < 25 per cent in the 2.8 < age ^4.6 Gyr age bin. We have assigned 39 N04 stars in this age bin to the Hercules stream (filled 
circles) using the Hercules stream UV phase-space (defined by F05, see Fig. [3} and the Hercules stream W phase-space (defined 
by F05, see Figs. [4] and [5jl. 
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Figure 2. U (left), V (middle) and W (right) velocity distributions of the N04 single stars with relative age errors < 25 per cent 
in the 2.8 < age 4.6 Gyr age bin (thin line) and a subsample of these stars that excludes the Hercules stream (thick line), as 
defined in Figs. [3] [4] and [5] 



3.2 



Comparison of N04 and F05 kinematic 
analysis 



F05 also published a VizieR catalogue (|Famaev et al.ll2004l ) 
accompanying their paper, based on the CORAVEL radial 
velocities of 6,691 northern hemisphere stars. Unfortunately, 
unlike the N04 stars, the F05 stars do not have individual 
stellar ages so F05 cannot be used to re- visit the age- velocity 
dispersion relation of the smooth background and investi- 
gate the effects of each dynamical stream on the relation. 
Instead in this section, we describe how F05 decomposed 
their sample into its constituent kinematic components and 
compare their phase-space substructure with N04 in order 
to physically motivate our assessment of Figs. [T] and [2] (and 
[XT] and Eg) , 

The basic technique to calculate space velocities is to in- 
vert the parallax to estimate the distance. This is then used 
to compute the tangential space motion components from 
the proper motions, which in combination with radial veloc- 
ities can be transformed to the Galactocentric frame. Due to 
the non-linearity of the parallax-distance transformation, if 



a-n/n is high, the inverse parallax is a biased estimator of the 
distance and an individual error on a space velocity cannot 
be derived from a simple first-order linear propagation of 
cr„. As all the stars in N04 are relatively nearby, most have 
parallaxes with ov/V better than 10 per cent and nearly all 
are better than 20 per cent. The availability of photometric 
distances with an uncertainty of only 13 per cent means the 
parallax inversion derived distances adopted for N04 stars 
when a^/n < 13 per ce nt only have a very small distance 
bias (|Brown et al.|[l997l ). which does not affect the calcula- 
tion of space velocity errors. 

The resulting space velocity diagrams, restricted to the 
8,589 single N04 stars, are plotted in the left columns of Figs. 
[31lland[l](c.f. N04's fig. 20 of the whole sample). While the 
U — W (Fig. [4} and V — W (Fig. E) diagrams show a smooth 
distribution, the U — V diagram (Fig. [3} shows abundant 
phase-space substructure. It consists of four tilted branches 
aligned along approximately constant V velocities that re- 
semble classic moving groups or stellar streams named, from 
top to bottom, the Sirius-UMa, Coma Berenices (or local), 
Hyades-Pleiades and £ Herculis branches. These branches 
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Figure 3. U - V diagrams of the 8,589 N04 single F-G dwarfs with space velocities (left column) and the 6,030 F05 K-M giants with 
space velocities, including spectroscopic binaries with centre-of-mass radial velocities (right column). The overdensity at (U,V) = 
(-45,-20) km s -1 in the left column is the Hyades open cluster (only in the N04 sample). In the bottom left panel, the N04 stars are 
colour-coded according to our assignment of N04 stars to the Hercules stream (green) and high-velocity stars (blue) in the Hercules 
stream UV phase-space (defined by F05, bottom right panel) but outside the Hercules stream W phase-space (see Figs. l4l and [5t. 
In the bottom right panel, the F05 stars are colour-coded according to their maximum likelihood base group assignment to known 
kinematic features of the solar neighbourhood: smooth background (black), Sirius moving group (magenta), young kinematics 
(yellow), Hyades-Pleiades supercluster (red), Hercules stream (green) and high- velocity stars (blue). 



are only loosely defined as overdensities in phase-space. 
N04's straight-forward space velocity derivation does not in- 
clude assignment of individual stars to these branches, only 
its position in phase-space. 

F05 binaries for which no centre-of-mass radial velocity 
could be estimated were removed from their sample using 
their binarity flag (B). A star was excluded in their kine- 
matic analysis if it was a spectroscopic binary (B = 0) or a 
visual binary with no orbit available (B = 5) or an uncer- 
tain case, either a spectroscopic binary or supergiant (B = 
3) or for a visual binary (B = 8), leaving 6,030 stars (5,311 
K and 719 M giants). Because the F05 giants are at dis- 
tances approximately an order of magnitude greater than 
the N04 dwarfs, the F05 distribution of a-n/n is correspond- 
ingly less accurate than that of N04. There are only 786 F05 
stars with a-n/n < 10 per cent, which is too small a sample 
to characterize the K and M giants in the solar neighbour- 
hood. ov/V < 20 per cent increases the sample size to 2,774 
stars (2,524 K and 250 M giants). With this accuracy, the 



distance bias is still very small and the first-order approx- 
imation for the calculation of space velocity errors is still 
valid. The resulting density of stars in the U — V diagram 
with ov/V < 20 per cent reveals overdensities in the same 
region of phase-space as the N04 branches (see F05's fig. 7). 
We therefore consider that the F05 kinematic phase-space 
structure seen at distances of a few hundred pc is consistent 
with that seen locally by N04. 

To make full use of all their 6,030 stars and to assign 
individual stars to each phase-space branch in order to in- 
vestigate their kinematics, F05 applied the Luri-Mennessier 
method to their data. This maximum-likelihood method is 
based on a Baye sian approach and is described in detail in 
iLuri et all (|l996l ). It requires a model describing the basic 
morphological characteristics of the sample (spatial distri- 
bution, kinematics and My) and a model of the selection 
criteria used to define the sample. A total distribution func- 
tion is built using these models to describe the observational 
characteristics of the sample. This consists of a linear com- 
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Figure A.U-W diagrams of the 8,589 N04 single F-G dwarfs with space velocities (left column) and the 6,030 F05 K-M giants with 
space velocities, including spectroscopic binaries with centre-of-mass radial velocities (right column). The overdensity at (U,W) = 
(-45,0) km s — 1 in the left column is the Hyades open cluster (only in the N04 sample). In the bottom left panel, the N04 stars 
are colour-coded according to our assignment of N04 stars to the Hercules stream (green) and high-velocity stars (blue) in the 
Hercules stream UV phase-space (denned by F05, see Fig. [3]l but outside the Hercules stream W phase-space (defined by F05, 
bottom right panel). In the bottom right panel, the F05 stars are colour-coded according to their maximum likelihood base group 
assignment to known kinematic features of the solar neighbourhood: smooth background (black), Sirius moving group (magenta), 
young kinematics (yellow), Hyades-Pleiades supercluster (red), Hercules stream (green) and high- velocity stars (blue). 



bination of partial distribution functions, each of which de- 
scribes a group of stars called a 'base group'. Each of these 
combines a velocity ellipsoid with a Gaussian magnitude dis- 
tribution and exponential height distribution perpendicular 
to the Galactic plane. The advantage of this phenomeno- 
logical model is that it can identify and quantify different 
subgroups in the data, which cannot be easily parametrized 
using other methods. A maximum-likelihood fit of the model 
parameters to the sample yield the model parameters that 
best represent the sample given the a priori models assumed. 

The expected space velocities deduced by the Luri- 
Mennessier method of the 6,030 F05 stars are plotted in the 
right columns of Figs. EHU and E](c.f. F05's fig. 9). Modelling 
this sample with a single base group returns a low maximum 
likelihood, suggesting to F05 that a single Schwarzschild el- 
lipsoid does not fit the kinematic properties of giant stars in 
the solar neighbourhood. The first acceptable solution they 
found required three base groups: bright giants/supergiants 
with 'young' kinematics, high-velocity stars and 'normal' 



stars. This last group contained further small-scale struc- 
ture, which F05 successfully modelled by three more statis- 
tically significant base groups. Although adding more groups 
produced statistically better solutions, F05 found that these 
solutions were not stable because they depended too much 
on the observed values of some individual stars. F05 identi- 
fied the resulting six base groups with known kinematic fea- 
tures in the solar neighbourhood (colour-coded in the bot- 
tom right panels of Figs. [3] [4] and [5}, which are also seen in 
the N04 phase-space. 

Despite all the differences described above, Figs. [3] 13] 
and [5] show the N04 and F05 phase-space is remarkably sim- 
ilar, suggesting the F05 kinematic base group assignment 
describes structures covering several hundred pc in the solar 
neighbourhood at least to first order, again suggesting the 
structures are not remnants of single small star clusters. We 
therefore consider the N04 branches to have the same physi- 
cal origin as postulated by F05 for their corresponding base 
groups. 
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Figure 5. V - W diagrams of the 8,589 N04 single F-G dwarfs with space velocities (left column) and the 6,030 F05 K-M giants with 
space velocities, including spectroscopic binaries with centre-of-mass radial velocities (right column). The overdensity at (V\W) = 
(-20,0) km s —1 in the left column is the Hyades open cluster (only in the N04 sample). In the bottom left panel, the N04 stars 
are colour-coded according to our assignment of N04 stars to the Hercules stream (green) and high-velocity stars (blue) in the 
Hercules stream UV phase-space (defined by F05, see Fig. [3]l but outside the Hercules stream W phase-space (defined by F05, 
bottom right panel). In the bottom right panel, the F05 stars are colour-coded according to their maximum likelihood base group 
assignment to known kinematic features of the solar neighbourhood: smooth background (black), Sirius moving group (magenta), 
young kinematics (yellow), Hyades-Pleiades supercluster (red), Hercules stream (green) and high- velocity stars (blue). 



3.3 Comparison of N04 and F05 phase-space 

We compare the space velocity distributions of the N04 and 
F05 stars in Fig. [6] This highlights the departures from 
Gaussianity of the U and V distribution functions caused 
by the presence of streams compared to W, which approx- 
imates a Gaussian, due to the tilt of these streams in the 
[/ V-plane in Fig. [3] me asured to be V = constant — 0A7U 
bv lSkulian et al.l lfl999l ). A direct comparison between the 
N04 and F05 distributions is interesting because Figs. [3] 
[4] and [5] show that the N04 substructures are more promi- 
nent than the F05 substructures compared to their smooth 
backgrounds apparently due to the difference in scaleheight 
between the two samples, illustrated in Fig. [7] The much 
lower scaleheight of the N04 sample and its more prominent 
substructure means dynamical streams are more prominent 
in stars close to the plane due to a mix of dynamics and 
sample age differences. 

Another difference is that there are similar numbers of 
stars in each of the F05 streams but the N04 V distribution 



in Fig.[6]shows that this is not the case in the N04 sample as 
it is numerically dominated by the much narrower peak of 
the Hyades-Pleiades stream at V ~ —20 km s _1 than in the 
F05 V distribution. This is because the N04 completeness 
out to ~ 40 pc in cludes the Hyades clu ster itself (visible in 
Figs. HH and 0. iFamaev et al. | l|2007l ) showed that in ad- 
dition to bound Hyades cluster stars, the Hyades-Pleiades 
stream is mainly composed of field-like stars but also partly 
of coeval stars evaporated from the Hyades cluster. F05 ac- 
tually propose that because the cluster is young (600 Myr), 
it should not have crossed many spiral structures, suggesting 
that the same transient spiral wave that formed it (by boost- 
ing star formation in its primordial gas cloud) could at the 
same time have given it and nearby field stars their peculiar 
velocity in the UV-pl&ne. For any evaporated cluster stars 
to get into the much higher scaleheight of the F05 sample 
requires them to have experienced phase- mixing and vertical 
disc heating. They are more likely to have had their initial 
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Figure 6. U (left), V (middle) and W (right) velocity distributions of the 8,589 N04 single F-G dwarf stars with space velocities 
(top row, thin black lines) and the 6,030 F05 K-M giants with space velocities (bottom row, thin black lines), including spectroscopic 
binaries with centre-of-mass radial velocities, colour-coded according to their maximum likelihood base group assignment to known 
kinematic features of the solar neighbourhood: smooth background (thick black lines), Sirius moving group (magenta), young 
kinematics (yellow), Hyades-Pleiades supcrcluster (red), Hercules stream (green) and high- velocity stars (blue). 



birth velocities erased and so are less likely to be identifiable 
as evaporated cluster stars in the F05 sample. 

F05 showed that their dynamical streams span several 
Gyr in age, which is supported theoretically by the simula- 
tions of |Pe Simone et al.l l|2004l ). Therefore, we expect the 
F05 streams in Fig. [6] to be present in the N04 sample over 
many age bins as narrow peaks in the V histograms in Figs. 
[2] and IA21 which allows us to trace the changing relative 
strengths of these streams with stellar age. This suggests 
that the majority of the peaks in the V histograms in Figs. 
[2] and IA2I are in fact dynamical streams rather than noise. 
We can also use Fig. |S] to interpret the N04 U and W dis- 
tributions. It suggests that the many of the minor peaks 
in these distributions are noise. However, like in V, the F05 
streams do not share similar mean U velocities so the N04 U 
distribution is also complex and non-Gaussian. We conclude 
that N04's treatment of the U and V distributions with a 
straight-forward dispersion was an inadequate parametriza- 
tion of these complex distribution functions. The interpre- 
tation of an age-velocity dispersion relation for U and V is 
non-trivial and beyond the scope of this paper. 

Nevertheless, these complex distribution functions raise 
the interesting issue as to what constitutes disc heating. In 



the introduction, we defined it as processes that change stel- 
lar velocities in any direction such that the velocity disper- 
sion increases. If this heating is random in direction it will fill 
and shape the velocity ellipsoid slowly over time. However, if 
processes like those responsible for dynamical streams (tran- 
sient spiral waves and/or resonances) overpopulate a specific 
subregion of phase-space, then the concentration of these 
stars relative to the background can actually change the 
overall components of the velocity dispersion by weighting 
them to the mean velocities of the streams. This scenario is 
illustrated in Table 1. It reveals that the F05 background 
dispersions are systematically slightly larger than the dis- 
persions of the background and thin disc streams combined. 
This shows that the F05 streams decrease the dispersions of 
the background and thin disc streams combined by weight- 
ing them to the mean velocities of the streams so that the 
combined dispersion differs slightly from the background dis- 
persion. Therefore, the mechanism responsible for the F05 
thin disc streams could be considered a disc 'cooling' agent. 
The inclusion of the Hercules stream and high- velocity stars 
increases the overall dispersion to greater than the back- 
ground alone so they could be considered as 'heating' agents. 
However, this classification is specific to the F05 sample. The 
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Figure 7. Three dimensional Galactic sky distribution of the 
8,589 N04 single F-G dwarfs with space velocities (grey) and 
the 6,030 F05 K-M giants with space velocities, including spec- 
troscopic binaries with centre-of-mass radial velocities, colour- 
coded according to their maximum likelihood base group as- 
signment to known kinematic features of the solar neighbour- 
hood: smooth background (black), Sirius moving group (ma- 
genta), young kinematics (yellow), Hyades-Pleiades supercluster 
(red), Hercules stream (green) and high- velocity stars (blue). The 
spherical polar axes have a radius of 900 pc, centred on the Sun. 
The viewer is in the Galactic plane, inside the solar circle at 
I = 30° (the Galactic centre is over the viewer's left shoulder). 
The viewer's angle has been chosen to illustrate the 63° between 
the Earth's equator and the Galactic plane (longest horizontal 
line), emphasizing that the F05 sample only contains stars vis- 
ible from the Earth's northern hemisphere. The nearest line of 
constant Galactic longitude is I = 30° (with the majority of 
the stars behind it). I increases to the right (anticlockwise) with 
lines of constant Galactic longitude every 30°. All the lines of 
constant Galactic longitude converge at the North Galactic Pole 
(6 = +90°, top) and at the South Galactic Pole (6 = -90°, bot- 
tom), with lines of constant Galactic latitude every 30°. 



Table 1. Comparison of N04 and F05 velocity dispersions (in km 
s -1 ) of combinations of different F05 kinematic groups: smooth 
background (B), Sirius moving group (S), young kinematics (Y), 
Hyades-Pleiades supercluster (HyPl), Hercules stream (He) and 
high-velocity stars. (Each dispersion is an outlier-resistant mea- 
sure of the dispersion about the centre of the velocity distribution 
calculated by using the median absolute deviation as the initial 
estimate, then each star is weighted using Tukcy's Biweight. For a 
Gaussian distribution, this is identical to the standard deviation.) 



Dispersion 




ay 


aw 


N04 all 


32.6 ± 0.2 


20.0 ± 0.2 


15.1 ± 0.1 


B 


33.8 ± 0.4 


17.4 ± 0.2 


17.8 ± 0.2 


Si, Y, HyPl 


22.3 ± 0.5 


12.2 ± 0.3 


7.0 ± 0.2 


B, Si, Y, HyPl 


31.8 ± 0.3 


16.4 ± 0.2 


15.7 ± 0.2 


B, Si, Y, HyPl, He 


33.7 ± 0.3 


19.6 ± 0.2 


15.8 ± 0.1 


F05 all 


35.6 ± 0.3 


22.0 ± 0.2 


17.1 ± 0.2 



bins (see Fig. IA2|) . representing the old thin disc, are more 
Gaussian and thus can be considered relaxed. We conclude 
that N04's treatment of the W distribution with a straight- 
forward dispersion may not be an adequate parametrization 
of the young thin disc but it is adequate for the old thin 
disc. This does not affect their general result because aw 
must increase with time in order for the older thin disc to 
have a higher scaleheight than the young thin disc. It does 
not affect our aims here because we are only concerned with 
the functional form of the relation for the old thin disc. 

It has long been known that traditional moving groups 
are indistinguishable in the vertical direction because the 
W velo city distributio n function appears to be well phase- 
mixed (|Dehnenl I1998T ). This is theoretically expected as 
phase-mixing of vertical space motions is more efficient than 
for the horizontal space motions because, for a realistic disc 
profile, the vertical frequency is a strong function of vertical 
energy and the vertical dynamical time is a factor of two 
shorter than the horizontal orbital period. Consequently, W 
motion is decoupled from U and V motions, suggesting that 
dynamical streams should not affect a star's W velocity. Ta- 
ble 2 in F05 suggest that all the streams share similar mean 
W ({W}). We test this in the following section. 



dominance of the N04 streams over the N04 background 
complicates the U and V distributions more than in F05. 
The differing strengths with stellar age of the streams, makes 
it difficult to estimate whether the overall ajj and av will 
be heated or cooled by the presence of streams compared to 
the background as a function of time, without prior knowl- 
edge of the properties of the N04 background, which are 
non-trivial to derive given the dominance of the streams. 
Therefore, we are unable to generalize which specific mech- 
anisms constitute disc heating or cooling because it depends 
on the sample considered. 

N04 excluded the youngest age bin from their power 
law fit to avoid biases due to unrelaxed young structures. 
The non-Gaussianity of some of the W histograms of the 
young thin disc (0.5 < age ^ 2.8 Gyr) bins (see Fig. IA2|) 
suggests that the N04 sample does not vertically relax until 
the stars are older than this age range. The W histograms of 
the 2.8 < age sC 4.6 Gyr (see Fig.[2j and 4.6 < age s: 7.6 Gyr 



3.4 Stellar warp of the Galactic disc 

iDehnenl l|l998T ) found that his outermost contours of phase- 
space density are skewed relative to the W axis in the sense 
that more stars moving faster than the Local Standard of 
Rest (LSR, +V) are also moving upward with respect to 
the LSR (+W) than downwards (— W). His fig. 6 shows 
that (W) is approximately constant at the VFlsr value for 
V < 10 km s^ 1 but curves upward for V > 10 km s . He 
interprets this curve as the signature of a stellar warp in the 
Galactic disc, which starts at the solar circle. If there are 
lots of N04 stars in the warp and these stars are included 
in our analysis, the departure from a constant (W) would 
artificially increase aw- This would bias the true aw so that 
it could appear as though there is more disc heating than 
in the assumed flat disc case (approximately constant VFlsr 
with respect to V). 

To check whether this source of aw bi as is pre s ent in 
the N04 data set, we reproduce fig. 6 from IDehnenl (|l998l ) 
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Figure 8. Robust mean vertical m otion ((W)) a s a function of V (left column, same V plot range, velocity bin centres and 2.5 km 
s — 1 velocity bin widths as fig. 6 in [Dolmen 1993) an d as a function of Galactic longitude, / (right column) for the N04 single F-G 
dwarfs and the F05 K-M giants. 8,386 N04 stars are included in the top left plot out of a possible 8,589 (98 per cent), whereas all 
are included in the top right plot. 5,852 F05 stars are included in the bottom left plot out of a possible 6,030 (97 per cent), whereas 
all are included in the bottom right plot. The triangle indicates the velocity of the Local Standard of Rest (LSR) with respect to 
the Sun: (VlSR,Wlsr) = (-Vq-Wq) = (-5.25,-7.17) km s" 1 jDehnen fc Binnevl Il998h . The horizontal solid lines along W L SR 
out to the LSR represents a flat Galactic disc out to the solar circle. The horizontal dashed lines outwards fro m the LSR rep resent 
a model of the flat disc continuing outside the solar circle instead of a stellar warp in the disc, as proposed bv lDehnenl d 19981 ). The 
X 2 tests compare the data to the model to the right of the vertical dotted lines at V = 15 km s — 1 (set to exclude the maximum V 
bins containing Sirius dynamical stream stars). Data points without error bars on (W) only contain one star and are thus excluded 
from the \ 2 tests (as is the empty F05 bin 42.5 < V ^45.0 km s _1 ). 



for the N04 and F05 samples in the left column of Fig. [8] 
Contrary to his claim, we do not find a signature of a stellar 
warp in the Galactic disc in either data set. Fig. shows 
that the number of stars decreases with +V in both data 
sets, which is reflected in the large error bars in Fig. [8] at 
high values of +V. 

The N04 sample in the top left panel of Fig. [8] is av- 
eraged over the entire sky. HI observations show that the 
Galactic disc is fiat out to approximately the solar circle 
and then the warp turns up towards the North Galactic 
Pole in the I ~ 90° direction, while turning towards the 
South Galactic Pole in the I ~ 270° direction. If the N04 
and F05 stars share the gas kinematics of the warp, then 
their participation in the warp would be seen by (W) > 
towards / ~ 90° and (W) < towards / ~ 270°. It is con- 
ceivable that if the warp is symmetric in amplitude about 
the plane and symmetric in I, then its signature could cancel 
itself out and result in a flat (W) as a function of V in the 
top left panel of Fig. [8] Plotting (W) as a function of / in 
the right column of Fig. [8] confirms this is not the case and 



Table 2. x 2 goodness-of-fit test results between the data and 
model in Fig. [8] n is the number of data points included in the 
model fit. n c is the number of constraints on each model fit (only 
one to normalize each model to each data set). The number of 
degrees of freedom is v = n — n c . The x 2 probability distribu- 
tion function for v degrees of freedom, Q(x 2 | I/ )j is t ne standard 
statistical significance of the x 2 test (Q(x 2 \ l/ ) = 1 — P(x 2 \v)). 



Sample 


n 


n c 


V 


x 2 


X> 


Q(xV) 


N04 V > 15 


13 


1 


12 


19.0 


1.58 


0.089 


F05 V > 15 


11 


1 


10 


16.5 


1.65 


0.087 


N04 all I 


36 


1 


35 


45.2 


1.29 


0.117 


F05 all / 


36 


1 


35 


47.1 


1.35 


0.083 



that there is no sinusoidal signature of the warp in I. The 
F05 sample becomes noisy at I > 230° because of incomplete 
coverage of the Galactic sky (see Fig. [7] and F05 fig. 6) . 

Table 2 shows that Q(xV) > 0.003 (-P(xV) < 0.997) 
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for all the tests so both data sets accept the flat disc model 
at the 3cr statistical significance level. However, the model is 
not strongly accepted as would be indicated by high Q(x 2 \ u ) 
values. This is because for each data set, there is one V bin at 
similar V values (32.5 < V < 35.0 and 35.0 < V < 37.5 km 
s" 1 for F05 and N04 respectively) for which the (W) is more 
than 2A{W) away from the flat disc model of constant Wlsr 
(i.e. no warp). It is these bins that cause the low Q(x 2 \ u ) 
values in Table 2. Respectively, the 6 and 4 stars in these 
bins do not share similar W velocities as would be expected 
if they were part of a warp or stream. These bins have very 
wide ranges of W velocities, perhaps because one or more of 
these stars are members of the thick disc or halo. The small 
number of stars in each bin results in the data points being 
offset from the rest. 

Therefore, contrary to iDehnenl l| 19981 ). we do not find 
a signature of the stellar warp in the Galactic disc in the 
N04 or F05 samples. This means that the N04 sample has 
an approximately constant {W). The small (~ 2 km s _1 ) 
excursions from Wlsr as a function of V do coincide with the 
V values of the streams in the N04 sample. These small net 
vertical motions are presumably due to the young unrelaxed 
stars evaporating from star clusters rather than field stars 
caught in the dynamical streams. 

The lower scaleheight of the N04 sample in Fig.[7]means 
its smooth background aw will be le s s tha n the F05 back- 
ground aw in Table 2. iFamaev et ail (|2007h showed that ~ 
15 per cent of the N04 Hyades-Pleiades stream is evaporated 
coeval Hyades cluster stars. Assuming an initial velocity dis- 
persion (Jo ~ 10 km s _1 , the evaporated Hyades stars will 
manifest themselves as a clump in the velocity distribution 
function that is narrow in V (because all the stars have iden- 
tical angular momentum) but has size > ao in U and W (?). 
Therefore, the N04 Hyades-Pleiades stream should have a 
larger aw than the purely dynamical F05 Hyades-Pleiades 
stream (this may also be the case for the other streams). 
Hence, the N04 background aw is likely to be more simi- 
lar to the N04 thin disc streams aw than in the F05 case. 
Thus the overall N04 W distribution is a sum of similar aw 
and similar (W) Gaussians, which is a Gaussian with sim- 
ilar aw and similar (W). Having established that the only 
component of the age-velocity dispersion relation that can 
be simply measured locally is the vertical component, the 
rest of the paper concentrates on re-visiting N04's age-aw 
relation. 



3.5 Reproducing the N04 age-aw~ relation 

Fig. [6] shows that the Hercules stream has a larger aw than 
the other F05 streams. This may be due to its different dy- 
nam ical origin from the o ther streams: perhaps the Galactic 
bar. iBensbv et ail (|2007T ) showed that the Hercules stream 
includes both younger and old disc stars and they suggest 
some thick disc stars may be included. The Hercules stream 
aw may be larger due to the thick disc stars in the stream 
retaining their characteristically higher aw , which is a prop- 
erty of the thick disc. Alternatively, the higher aw may be 
a dynamical signature of the bar, imparted to both its con- 
stituent thin and thick disc stars as they interacted with the 
bar. If it is the latter and the N04 Hercules stream aw is 
sufficiently different from the N04 smooth background and 
thin disc streams, then the thin disc Hercules stream stars 



could bias the age-ay^ relation. The smooth W distributions 
in Figs. [4] and [5] suggest aw is a valid description of the W 
distribution function, unbiased by kinematic substructure 
because the W distribution is well phase-mixed. By deriv- 
ing the age-ay^ relation both with and without the Hercules 
stream, we can simultaneously test whether the relation is 
biased by the stream and whether the W distribution is well 
phase-mixed, implying the presence of the thin disc streams 
does not affect the relation. 

We use the most extreme F05 Hercules stream stars in 
U and V to define the N04 Hercules stream UV region in 
Fig. 13] Figs. |3] and [S] are also used to identify high- velocity 
stars within the UV region of the Hercules stream that have 
W velocities that are outside the Hercules stream W region 
and thus are more likely to be part of the thick disc velocity 
ellipsoid. This technique will only identify high- velocity stars 
with more extreme thick disc or halo kinematics. Therefore, 
the removal of Hercules stream stars will also remove some 
high- velocity stars but this is of no consequence to our study. 
These stars are colour-coded in Figs. [3] [4] and [5] Every fig- 
ure from here onwards both includes and excludes Hercules 
stream stars. 

Fig. [9] shows that a robust (bisquare weighted) mean 
age of an age bin is sometimes quite different to the central 
age of that bin. By comparison with N04's fig. 31, we find we 
are best able to reproduce their age-er^ relation using both 
robust ages and robust dispersions in Fig. [TUJ The error bars 
in Fig. [TO] are the errors in calculating a standard deviation: 



Aaw 



aw 



(1) 



where N is the number of stars in each age bin used to 
calculate aw- The logarithmic spacing of the age bins in 
Fig. [TU] is a natural consequence of N04's decision to have 
equal numbers of stars in each of their age bins because the 
distribution of stellar ages is logarithmic. This also naturally 
leads N04 (and us) to fit a power law: 



aw 



a age 



(2) 



where a is a constant and k is the scaling (heating) exponent. 
Taking the log of both sides of equation [5] gives 



logio(crw) = fclog 10 (age) + log 10 (a), 



(3) 



which produces a linear relationship where k is now the slope 
and logio(a) is the intercept. A linear least-squares fit gives 
k — 0.48 ± 0.26 (including the Hercules stream, k = 0.50 ± 
0.25 excluding Hercules stream), which is similar to N04's 
k — 0.47±0.05. However, we can only reproduce their quoted 
uncertainty by excluding errors from the fit: k — 0.45±0.04. 
We are thus able to reproduce and verify the N04 age-aw 
relation using the rounded up catalogue data. 

These values of k are close to the k — 0.50 predicted by 
lHanninen fc Flvnnl l|2002l ) for vertical disc heating by halo 
black holes of mass 1 x 10 7 Mq . The goodness of fit of these 
power laws are therefore a test of how well the data are rep - 
resented by this source of heating. [Hanninen fc Flvnnl (2002) 
also predict k = 0.26 for the vertical heating exponent by 
GMCs, where the constant a in equation[2]needs to be set by 
observations. We have plotted their power law prediction in 
Fig. 1101 setting a = ao — 10 km s _1 so that at the youngest 
age included in the fits (1.25 Gyr), aw ~ 10.6 km s _1 . 

The linear a,ge-aw plane in Fig. [11] highlights the un- 
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Figure 9. Discrete age distribution of the N04 single stars with relative age errors < 25 per cent, split between the youngest 7 age 
bins from Fig. IA2I (left) and the oldest 3 age bins from Fig. IA2I (right). The thin lines correspond to the whole sample and the thick 
lines are the subsample excluding the Hercules stream. Because the N04 catalogue stellar ages are provided to the nearest 0.1 Gyr, 
each 0.1 Gyr histogram bin only contains stars with the exact age that its histogram step is centred on e.g. the youngest histogram 
bin, centred on 0.6 Gyr, only contains stars with exactly 0.6 Gyr ages but the histogram step spans 0.55 to 0.65 Gyr. The vertical 
solid lines (including the vertical edges of the plots) are the boundaries to the age bins from Fig. IA2I The dotted lines are the 
central ages of each age bin. The robust mean ages within each age bin for the whole sample (thin dashed lines) and subsample 
(thick dashed lines) are indisguishable apart from the oldest age bin. 
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Figure 10. N04 logarithmic age-cr^v relation: robust aw of N04 single stars with relative age errors < 25 per cent in each of the 
10 age bins in Fig. IA2l as a function of the robust mean age of each bin including (left, filled triangles) and excluding (right, open 
triangles) the Hercules stream. As in N04, the youngest and oldest age bins have been excluded from the power law fit to avoid 
biases due to unrelaxed young structures and thick disc stars respectively. 



equal spacing in age of the age bins as a result of the N04 
procedure requiring the same number of stars in each bin. 
The densely distributed age bins at the young end confirm 
the presence of an age-<7w relation for the young thin disc. 
However, N04 do not consider whether any other models fit 
their data better than a power law. They assume a power 
law is the best- fitting model and interpret this as evidence of 
continued disc heating throughout its lifetime. The relative 
age error of < 25 per cent translates to sparse sampling of 
the old thin disc (two data points) and the thick disc (one 
data point). These are the only three age bins between 2.9 
and 15.2 Gyr and only the former two bins are included in 



the power law fit. Therefore, the power law is not as tightly 
constrained at the old end as it is at the young end. 

We test N04's assumption that a power law is the best- 
fitting model by generating alternative time dependencies 
of disc heating models (overlaid in Fig. 1 1 1 p and comparing 
them to the data. Our main aim is to establish whether the 
age-trw relation saturates for older stars. This means we 
are not so concerned with the exact functional form of the 
relation younger than a hypothetical saturation age. Thus 
we simply define models as linear fits of the data points up 
to a hypothetical saturation age, after which aw is constant 
with time, set by the aw at the hypothetical saturation age. 
The number of constraints on each power law model fit (n c ) 




Figure 11. N04 linear a,ge-aw relation: robust aw of N04 single stars with relative age errors < 25 per cent in the 8 age bins 
included in the power law fit in Fig. llOl as a function of the robust mean age of each bin including (left, filled triangles) and excluding 
(right, open triangles) the Hercules stream. The dashed lines are models representing no disc heating, where the constant saturation 
level is set by the dispersion in the youngest bin in Fig. 1101 The solid lines are linear fits to each set of data points up to 2.5, 3.5, 
4.5, 5.5 and 6.5 Gyr. Thereafter these models saturate at constant dispersions. 



is 2 because the number of degrees of freedom (v) is reduced 
by the two constraints of the fit: the slope and intercept. 
This is also the case for each linear fit before saturation 
model but v is further reduced by one more constraint: the 
choice of saturation epoch. The \ 2 statistic for our purpose 
of comparing the observed distribution of aw to a model is 

where a w is the observed aw of the ith age bin, measured as 
a robust aw, and a w is the robust aw expected according 
to a model. 

As expected, Figs. [10] and [11] show that the largest ef- 
fect in removing the Hercules stream from the N04 binning 
sample is seen in the oldest thick disc age bin, which is 
not included in the a,ge-aw relation. Excluding the Hercules 
stream stars removes Hercules stream thin disc stars from 
the thick disc region of phase-space and so the thick disc 
aw increases slightly. Figs. [TU] and [TT] and Table 3 show 
that removing the Hercules stream has a negligible effect 
on the &ge-aw relation. This may suggest that there are in- 
sufficient numbers of Hercules stars in each age bin for the 
hypothetically higher Hercules aw to affect the relation. It 
also suggests that dynamical streams are well phase-mixed 
in W, supporting our argument that the age-ow relation 
can be constrained from the N04 data. 

Table 3 shows that the N04 age-binned data (includ- 
ing the Hercules stream) accept the power law and thus 
halo black holes are not ruled out as vertical disc heating 
sources. As expected the data strongly reject a no disc heat- 
ing model and linear fits before saturation at 2.5 and 3.5 Gyr 
models. The data also reject the GMC heating power law. 
However, the Q{x 2 \ v ) values for this model are very sensi- 
tive to the initial aw- As mentioned earlier, if the unrelaxed 
early age bins actually have larger values of aw, then the 



lHanninen fc Flvnnl (|2002T ) model could perhaps be accepted 
by the data. We confirm one of N04's key results that disc 
heating does not saturate at an early stage. The linear fit 
before saturation at 4.5 Gyr model is only marginally re- 
jected by the data including the Hercules stream, while the 
data excluding the Hercules stream accept this model. Mod- 
els with linear fits before saturation at 5.5 and 6.5 Gyr are 
accepted by the data including and excluding the Hercules 
stream. Hence, the N04 data do not require a power law and 
actually statistically prefer saturation after ~ 4.5 Gyr. 

3.6 Higher resolution age binning 

Higher resolution age binning is required to try to further 
rule out other disc heating models. Linear age binning, which 
needs bins of constant age range, will sample the old thin 
disc with more data points. Selecting the range of each age 
bin to be 1 Gyr wide more than doubles the number of data 
points sampling the old thin disc and thick disc. It also more 
than halves the number of age bins sampling the young thin 
disc but this regime has already been well constrained in the 
previous section. 

We also chose an absolute age error of < 1 Gyr so that 
any star can only shift laterally by a maximum of one age 
bin. This is not the case with the N04 binning. The prepon- 
derance of young stars in the N04 magnitude-limited sam- 
ple (because of their higher luminosity) means their bins 
are more narrow than the typical age errors (< 25 per cent) 
so the errors of some young stars span more than one age 
bin. Absolute age errors < 1 Gyr means selecting stars with 
°"ag| h - a g e < 1 Gyr and age - cr™ < 1 Gyr. Selecting sin- 
gle stars with absolute age errors < 1 Gyr actually samples 
more stars (3,083) than the relative age error of < 25 per 
cent (2,801) because the vast majority of N04 stars are be- 
tween 1 and 3 Gyr old (see N04's fig. 17). The youngest and 
oldest stars with an absolute age error < 1 Gyr are 0.6 and 
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Table 3. x 2 goodness-of-fit test results between the N04 age binning data and disc heating models in Fig. 1111 n is the number 
of data points included in the model fit. n c is the number of constraints on the model fit. The number of degrees of freedom is 
v = n — n c . The x 2 probability distribution function for u degrees of freedom, Q(x 2 \ v )> i s the standard statistical significance of the 
X 2 test (Q(x 2 \f) = 1 - P{x 2 W))- If Q(X 2 W) > 0-003 (P{x 2 \v) < 0.997), the data accept the model at the 3ct statistical significance 
level. If Q(x 2 \ u ) < 0.003 (P(x 2 \v) > 0.997), the data reject the model at the 3cr statistical significance level. 
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Figure 12. Discrete age distribution of the N04 single stars with absolute age errors < 1 Gyr, split between the youngest 3 age 
bins (left) and the oldest 7 age bins (right). The thin lines correspond to the whole sample and the thick lines are the subsamplc 
excluding the Hercules stream. Because the N04 catalogue stellar ages are provided to the nearest 0.1 Gyr, each 0.1 Gyr histogram 
bin only contains stars with the exact age that its histogram step is centred on e.g. the youngest histogram bin, centred on 0.6 Gyr, 
only contains stars with exactly 0.6 Gyr ages but the histogram step spans 0.55 to 0.65 Gyr. The vertical solid lines (including the 
vertical edges of the plots) are the boundaries to the age bins. The dotted lines are the central ages of each age bin. The robust 
mean ages within each age bin for the whole sample (thin dashed lines) and subsample (thick dashed lines) are indisguishable until 
the three oldest age bins. 



10.5 Gyr, constraining our age interval to 0.5 < age ^ 10.5 
Gyr, giving us ten 1 Gyr wide age bins. Fig. [T^] shows the 
age distribution of our new sample. 

Figs. [13] and [14] show the space velocity diagrams and 
histograms respectively of a representative age bin from our 
new sample (plots of all our age bins are in the Appendix) 
that is most similar in age to the N04 age bin in Figs. [T] 
and E] Unlike the N04 binning in Fig. (and Fig. QOJ), the 
age bin in Fig. [14] (and the two younger bins, all representing 
the young thin disc, in Fig. lA4Jl do not suffer from statistical 



noise. They show the non-Gaussianity of the U and V dis- 
tributions compared to W more clearly than in Figs. [5] and 
IA2I Again, the young thin disc W histograms do not look 
completely relaxed. Fig. IA4I clearly shows that only includ- 
ing stars with errors < 1 Gyr severely reduces the numbers 
of stars in the older bins. The oldest W distribution that 
is still approximately a complete symmetric Gaussian and 
statistically significant with 26 stars is in bin 5.5 < age ^ 
6.5 Gyr. Therefore, after the removal of statistically insignif- 
icant age bins older than 6.5 Gyr, our higher age resolution 
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Figure 13. U — V (left), U — W (middle) and V — W (right) space velocity diagrams of the 527 N04 single stars with absolute age 
errors < 1 Gyr in the 2.5 < age $C 3.5 Gyr age bin. We have assigned 62 N04 stars in this age bin to the Hercules stream (filled 
circles) using the Hercules stream UV phase-space (defined by F05, see Fig. |3j and the Hercules stream W phase-space (defined 
by F05, see Figs. [4] and [5J|. 



W 

O 




-80 -40 40 80 -80 -40 40 80 -80 -40 40 80 
U (km s" 1 ) V (km s" 1 ) W (km s" 1 ) 



Figure 14. U (left), V (middle) and W (right) velocity distributions of the N04 single stars with absolute age errors < 1 Gyr in 
the 2.5 < age sC 3.5 Gyr age bin (thin line) and a subsample of these stars that excludes the Hercules stream (thick line), as defined 
in Figs. [3] Hand [5] 



binning only adds one extra bin with which to constrain the 
age-aw relation of the old thin disc. The 5.5 < age ^ 6.5 
Gyr is the oldest age bin included in the power law fit in 
Fig. 1151 The power law fit to our binning produces a higher 
scaling (heating) exponent than the N04 binning power law: 
k = 0.63 ± 1.21 including the Hercules stream (k = 0.60 ± 
1.22 excluding the Hercules stream). 

We plot all the age bins from Fig. [12] in Fig. [15] but dif- 
ferentiate between the statistically significant and insignifi- 
cant bins. Although not included in the power law fit, the 
oldest age bin in N04's binning in Fig. [10] (and N04's fig. 31), 
seems to be a continuation of the power law, suggesting it is 
part of a continuing heating process. The oldest age bin in- 
cludes some thick disc stars, which have been shown by many 
authors to have approximately double the thin disc aw- This 
thick disc signature is not seen in the N04 binning because 
the oldest age bin spans the oldest part of the thin disc and 
all the thick disc so it is numerically dominated by old thin 
disc stars that decrease aw- Our higher age resolution bins 
better resolve the thin disc-thick disc transition region. The 
age bins centred on 8 and 9 Gyr in Fig. [15] show a mod- 



est increase in aw compared to the power law fits, reaching 
values more associated with the thick disc, 35 ± 3 km s~ 



|Chiba fc BeersfeOOOT ). 39±4 km s" 1 l|Soubiran et al.ll2003l ). 
than the thin disc. This is in agreement with the jumps in 
the age- veloci ty dispersion relation at 8 and 9 Gyr fo u nd re - 
spec tively by Sornmer-Larsen fc Antonuccio-Delogrj (| 19931 ) 



and iQuiilen fc Garnetd ( 200 ll ). the latter being significant 
at the 2-3cr level. The abrupt increase is interpreted as a 
formation signature of the thick disc, which is consistent 
with models where the thick disc was formed by dramatic 
disc heating when a satellite galaxy falls o nto the Galac- 
tic disc in a minor merger ~ 10 Gyr ag o ( ? l Freeman! Il99ll : 
iToth fc Ostrikerll 19921 : buinn et al.l(l993l ). Despite their sta- 
tistical insignificance, it is tempting to associate these bins 
with the formation of the thick disc. 



Again as expected, Figs.[l5]and[T6]show that the largest 
effect in removing the Hercules stream from our binning 
sample is seen in the oldest thick disc age bins, which are 
not included in the a,ge-aw relation. Figs. 1151 116l and Table 
4 show that removing the Hercules stream has a negligible 
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Figure 15. Our logarithmic age-criy relation: robust cryy of N04 single stars with absolute age errors < 1 Gyr in each of the 10 age 
bins in Fig. IA4l as a function of the robust mean age of each bin including (left, filled triangles) and excluding (right, open triangles) 
the Hercules stream. As in N04, the youngest age bin has been excluded from our power law fit (long dashed line) to avoid biases 
due to unrelaxed young structures. The oldest age bin included in our fit is centred on 6 Gyr. Older age bins are excluded because 
they contain insufficient numbers of stars (see Fig. IA4t for their oyy values to be statistically significant (dashed error bars). The 
distribution of W velocities in the age bins centred on 10 Gyr do not yield robust values of crw - Therefore, their standard deviation 
is plotted instead. 




Figure 16. Our linear age-cr^ relation: robust ayy of N04 single stars with absolute age errors < 1 Gyr in the 5 age bins included 
in our power law fit in Fig. ll5l as a function of the robust mean age of each bin including (left, filled triangles) and excluding (right, 
open triangles) the Hercules stream. The dashed lines are models representing no disc heating, where the constant saturation level 
is set by the dispersion in the youngest bin in Fig. 1151 The solid lines are linear fits to each set of data points up to 2.5, 3.5, 4.5, 
5.5 and 6.5 Gyr. Thereafter these models saturate at constant dispersions. 



effect on the age-aw relation, allowing the relation to be 
constrained by our age binning as well as the N04 binning. 

The same disc heating models as before are considered 
with the results in Table 4. Again, we have plotted the GMC 
power law in Figs. 1151 and [T6l setting a = ao = 10 km s _1 
so that at the youngest age included in the fits (1.55 Gyr), 
aw ~ 11.2 kms" 1 . With our binning there are fewer degrees 
of freedom, so the Q(x 2 \ l/ ) values are generally smaller in 
Table 4 than in Table 3 for the N04 binning. Our binned data 



including the Hercules stream accept our power law model 
but marginally reject it when the stream is excluded. Both 
including and excluding the Hercules stream in our binned 
data also marginally reject the N04 power law fit, strongly 
reject a no disc heating model, reject the GMC power law 
(although as discussed earlier, this may not be conclusive), 
reject a linear fit before saturation at 3.5 Gyr model but 
accept models with linear fits before saturation at 4.5, 5.5 
and 6.5 Gyr. Hence, our binning reinforces results from the 
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Table 4. \ 2 goodness-of-fit test results between our age binning data and disc heating models in Fig. 1161 n is the number of data 
points included in the model fit. n c is the number of constraints on the model fit. The number of degrees of freedom is v = n — n c . 
The x 2 probability distribution function for v degrees of freedom, Q(x 2 \ u )i is th e standard statistical significance of the x 2 test 
{Q(X 2 W) = 1 - P (X 2 W))- If Q(X 2 W) > 0.003 (P(x 2 W) < 0.997), the data accept the model at the 3cr statistical significance level. 
If Q{x 2 \ u ) < 0.003 (P(x 2 \ u ) > 0.997), the data reject the model at the 3<r statistical significance level. 
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previous section that the data do not require a power law 
and actually statistically prefer saturation at ~ 4.5 Gyr. 



4 DISCUSSION 

This paper re-visits the Galactic thin disc age-velocity dis- 
persion relation. We show that dynamical streams compli- 
cate the U and V velocity distribution functions such that a 
straight-forward dispersion is not an adequate parametriza- 
tion of their distribution functions. This calls into question 
the appropriateness of power law fits to the au and ay re- 
lations with stellar age. It shows that observationally mea- 
sured err/ and ov can only be used to constrain in-plane 
disc heating models if the N04 sample can be kinematically 
decomposed (like in F05) so that only the smooth veloc- 
ity ellipsoidal background is used, excluding the dynamical 
streams. Their strong presence in the solar neighbourhood 
may cast doubt on the applications of au in determining the 
asymmetric drift, the extrapolating to zero au to recover 
the V solar motion and the interpretation of the velocity 
ellipsoid dispersion axis ratios av/au, aw /au and aw lay. 
F05 have already pointed out using their sample that the 
underlying velocity ellipsoid, after the removal of streams 
is not centred on the commonly accepted radial anti-centre 
motion. Instead it is centred on (U) = —2.78 ± 1.07 km 
s - . Inclusion of the streams yields the commonly accepted 
value: (U) = -10.25 ± 0.15 km s" 1 . 

The &ge-aw relation can be constrained using the N04 
data set because the dynamical streams are well phase- 
mixed in W, as illustrated by removal of the Hercules stream 
having no effect on the relation, which also suggests that the 
Galactic bar does not vertically heat the local disc. We do 
not find any signature of the stellar warp in the Galactic 



disc so we do not need to exclude any N04 stars from the 
a,ge-aw relation on this basis, because they are all free from 
this potential source of aw bias. 

We reproduce the N04 a,ge-aw relation and confirm 
their finding that vertical disc heating does not saturate at 
an early stage (before ~ 4.5 Gyr). Their logarithmic age bin- 
ning naturally led them to fit a power law. However, they 
did not consider other disc heating models and so their con- 
clusion was that the power law fit was evidence that verti- 
cal disc heating is continuous throughout the lifetime of the 
thin disc. We have considered other disc heating models that 
saturate after different epochs. Our key new result is that a 
power law is not required by their data. A power law fit is 
statistically similar to disc heating models which saturate af- 
ter ~ 4.5 Gyr (including and excluding the Hercules stream). 
These saturation model fits are consi stent with the observa - 
tional turnover at ~ 5 Gyr found bv lCarlberg et aT, | (|1985| ). 
Simulat ions of vertical disc heati ng solely from GMCs find k 
= 0.26 (|Hanninen fc FlvnrJ|20o3 ). Figs.fTTIandfrolshow that 
this power law is consistent with a minimal increase in aw 
for old stars, not too dissimilar to the models that saturate 
after 3.5 and 4.5 Gyr. 

To see if the data could rule out any of the disc heating 
models considered, we examined the relation at a higher 
age resolution but found very similar results to the N04 age 
binning. Therefore, two diametrically different age binning 
methods using the best data available for constraining the 
&ge-aw relation still cannot distinguish between competing 
disc heating models. The issue is unlikely to be resolved until 
there are many more space velocities and ages measured for 
old thin disc star s. This may need to w ait until the ESA 
Gaia satellite era (|Perrvman et al.ll200ll ). 

Nevertheless, our results provide new constraints to be 
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fulfilled by disc heating simulations. Prior to this paper, 
new disc heating simulations would have tried to explain 
the continued disc heating found in all three cardinal di- 
rections by N04. Now, we have shown that new simula- 
tions cannot use in-plane age-velocity dispersion relations 
derived using stellar samples containing dynamical streams 
and straight-forward dispersions. Also, the new simulations 
must be able to reproduce the complex U and V sub- 
structure in the solar neighbourhood, caused by dynami- 
cal streams, while modelling in-pla ne disc heating mech- 
anisms (e.g. |Pe Simone et aLll2004T ). They must also take 
into account that these dynamical streams dominate the so- 
lar neighbourhood velocity distribution function and include 
evaporating star clusters, which may have been created by 
the spiral wave that spawned the streams. 

Assuming our c hoice of initial aw and the modelling of 
lHanninen fe Flvnnl (|2002h are both realistic, simulations of 
vertical disc heating solely from GMCs appear unable to re- 
produce the amplitude of our age-ow relation on their own. 
This suggests that another vertical disc heating mechanism 
is or has been involved in the evolution of the Galactic thin 
disc. We have shown that the data do not rule out heating by 
halo black holes. However, this option has no observational 
support. Black holes passing through the disc would proba- 
bly reveal themselves as high proper motion X-ray sources 
and no such objects have been found. 

In the context of the hierarchical galaxy formation 
paradigm, another vertical heating mechanism is minor 
mergers with dwarf galaxies. We find extremely tentative 
evidence of an abrupt feature in the age-trw relation at ~ 
8 Gyr, which could be the sudden vertical heating of an 
early thin disc int o the thick disc due to a major merger 
( CJuinn et alJll993h . The relation does not exhibit any other 
strong features like this, suggesting it is unlikely that the 
Milky Way suffered another major merger in the last ~ 8 
Gyr but it does not rule out minor mergers during this pe- 
riod of relatively quiescent disc evolution. In Figs. llll and ll6l 
the old thin disc bins with aw greater than the GMC power 
law may be the signature of mi nor mergers that too k place 
> 3 Gyr ago at redshift > 0.5. IVelazquez fe White! (|l999l ) 
showed that disc galaxies can accrete quite massive satel- 
lites without destroying the disc, particularly if the orbits 
are retrograde. They found that a satellite on a retrograde 
orbit with an initial mass of 20 per cent of the Milky Way's 
disc could increase the velocity ellipsoid at the solar neigh- 
bourhood by Aaw = 6 km s _1 . The same satellite, but on a 
prograde orbit, leads to Aaw = 12 km s" 1 . It is intriguing 
that the difference between the observed age-erw relation 
and the predicted relation only heated by GMCs is ~ 6 km 
s _1 at an epoch ~ 5 Gyr ago in Fig. 1161 The agreement 
between the observed age-aw relation and the predicted re- 
lation only heated by GMCs is good between the present 
day and ~ 3 Gyr ago, which may suggest that no minor 
mergers with significant satellite mass have occurred in the 
last ~ 3 Gyr. 

The RA dial Velocity Experiment (RAVE, 

ISteinmetz et al.l 12006( 1 is measuring the radial veloci- 
ties of stars further from the Sun than the N04 and F05 
samples. RAVE vertical velocities can be used to test the 
vertical smoothness of the velocity distribution further 
from the Galactic plane than possible with N04 and F05, 
as well as to search for tidal streams falling onto the Milky 



Way disc from satellite galaxies (Seabroke et al. 2007, in 
prep.), which may be responsible for some of the vertical 
disc heating observed in the age-ow relation in the solar 
neighbourhood. 
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APPENDIX A: SPACE VELOCITY 
DISTRIBUTIONS IN ALL THE AGE BINS 

Al N04 age bins 

For completeness, all the N04 age bins have their space ve- 
locity diagrams plotted in Fig. IA1I and their space velocity 
histograms plotted in Fig. IA2I 

A2 Our age bins 

For completeness, all our age bins have their space velocity 
diagrams plotted in Fig. IA3I and their space velocity his- 
tograms plotted in Fig. IA4I 

This paper has been typeset from a TjtX/ DTpX file prepared 
by the author. 
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Figure Al. U — V (left), U — W (middle) and V — W (right) space velocity diagrams of the N04 single stars with relative age 
errors < 25 per cent plotted in the 10 age bins from Fig, IA2I (youngest at the top and oldest at the bottom). We have assigned N04 
stars to the Hercules stream (filled circles) and high-velocity stars (open circles) in the Hercules stream UV phase-space (defined 
by F05, see Fig. [3} but outside the Hercules stream W phase-space (defined by F05, see Figs, l4l and [Bl, 
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Figure A2. U (left), V (middle) and W (right) velocity distributions of the N04 single stars with relative age errors < 25 per cent 
plotted in the 10 age bins that dissect their corresponding age distribution to give approximately equal numbers (N a [[) of stars in 
each age bin (thin lines) and a subsample (N gu u) of these stars that excludes the Hercules stream (as defined in Figs. [3] [4] and [5j 
plotted in the same age bins (thick lines). The boundaries of each age bin are given in Gyr. 



Re-visiting the relations: Galactic thin disc age-velocity dispersion relation 23 





-50 
-100 



-50 
-100 



-50 
-100 



-50 
-100 



-50 
-100 



-50 
-100 



-50 
-100 



-50 
-100 



-50 
-100 



-50 
-100 





50 









-50 




50 









-50 




50 









-50 




50 









-50 




50 







1 


-50 


in 




E 


50 











>■ 






-50 




50 









-50 




50 









-50 




50 









-50 




50 









-50 



O 





o 










> ' . ■ * ' 


s 


■ • ■ •*,*/' 

© 






■ 


■ 


• 













100 
U (km s" 1 ) 



100 -100 -50 
U (km s" 1 ) V (km s" 1 ) 



Figure A3. U — V (left), U — W (middle) and V — W (right) space velocity diagrams of the N04 single stars with absolute age 
errors < 1 Gyr plotted in the 10 age bins from Fig. 1141 (youngest at the top and oldest at the bottom). We have assigned N04 stars 
to the Hercules stream (filled circles) and high-velocity stars (open circles) in the Hercules stream UV phase-space (defined by F05, 
see Fig. [3} but outside the Hercules stream W phase-space (defined by F05, see Figs. l4l and [5t. 
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Figure A4. U (loft), V (middle) and W (right) velocity distributions of the N04 single stars with absolute age errors < 1 Gyr 
plotted in the 10 age bins that dissect their corresponding age distribution to give equal width (1 Gyr) age bins. The boundaries of 
each age bin are given in Gyr. N a jj is the total number of stars in each age bin (thin lines) and N—,1. is the number of these stars 
in the subsample that excludes the Hercules stream (as defined in Figs. [3] [4] and [5} plotted in the same age bins (thick lines). 



